Universal temperature dependence of the magnetization of gapped spin chains 
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Temperature dependence of the magnetization of the Haldane spin chain at finite magnetic field 
is analyzed systematically. Quantum Monte Carlo data indicates a clear minimum of magnetization 
as a function of temperature in the gapless regime. On the basis of the Tomonaga-Luttinger liquid 
theory, we argue that this minimum is rather universal and can be observed for general axially 
symmetric quasi-one-dimensional spin systems. Our argument is confirmed by the magnetic-field 
dependence of the spin-wave velocity obtained numerically. One can estimate a magnitude of the 
gap of any such systems by fitting the experimental data with the magnetization minimum. 

PACS numbers: 75.10.Jm 

I. INTRODUCTION 

The discovery of the Haldane gap 1 established an important paradigm followed by the extensive studies on gapped 
spin chains. In particular, recent remarkable progress in high magnetic field experiments provide various data on 
the closing of the gap by the applied magnetic field. They stimulate a renewed interest in corresponding theoretical 
studies. However, because of methodological reasons the theories on gapped spin systems still lie far behind those of 
the exactly solvable 5 = 1/2 Heisenberg chain 2 . Actually at present the only clue to understand the experiments are 
the phenomenological effective theories 3 -! 4 - and the numerical calculations which support them. 

In this paper, we study the temperature dependence of the magnetization of the S = 1 Haldane chain for wide 
range of magnetic field. We particularly focus on the cases where the magnetic field, h, exceeds the gap, A. The 
magnetization is one of the most fundamental quantity observed in experiments. In spite of the simplicity of the prob- 
lem, it is not easy to give a systematic understanding of the magnetization in such systems. Numerical approaches 
to the problem also had difficulties. For example, quantum Monte Carlo (QMC) method**^ had difficulty due to 
small acceptance ratio in a finite magnetic field. This was overcome very recent ly£ & 1 1 1 1 1 1 2 with new formulations 
such as Stochastic Series Expansion (SSE). Taking advantage of such developments, we revisit the problem by com- 
bining various analytical approaches with numerical results obtained with modern techniques (QMC with SSE, and 
density matrix renormalization group(DMRG)). We clarify universal features in the temperature dependence of the 
magnetization in the Haldane chain under an applied magnetic field. Our findings are applicable to general gapped 
one-dimensional spin systems, as well. 

The Hamiltonian of the Heisenberg spin chain is given by 

N N 

H = jJ2s j .S j+1 -h^2S], (1) 

where J > and N is the system size. The ground state of 5 = 1 Haldane chain is non-magnetic (a singlet state) and 
has a finite energy gap, A = 0.4105 J, which is estimated numericall y 13 ' 14 . The system undergoes a quantum phase 
transition at finite magnetic field, h = h c (= A); in the Haldane gap phase h < h c , the magnetization exponentially 
vanishes toward zero temperature, which is confirmed experimentally, e.g. in Ref . [TEl. On the other hand, in the h > h c 
phase (gapless regime) the magnetization is finite even at T = 0. This regime can be described by the Tomonaga- 
Luttinger (TL) liquid 3 - as in the 5 = 1/2 Heisenberg chain. Therefore, one might expect that its magnetization also 
follows that of the 5 = 1/2 Heisenberg chain, which is well-known as the Bonner-Fisher curved. This is, however, 
not the case. We show that the magnetization curve as a function of T in the 5 = 1 gapless regime has a minimum 
at T = T m and that the value of T m depends on the magnetic field. Although this behavior resembles that of the 
three-dimensional Bose-Einstein condensation (BEC) of magnons 1 ^, the origin obviously differs between the two cases 
since the BEC at T ^ never occur in the one-dimensional systems 1 ^. We prove that this phenomena can be observed 
in general one-dimensional spin systems (chains and ladders) with the rotational symmetry about the magnetic field 
direction (axial symmetry). 
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The magnetization of the S — 1 gapless phase has already been studied in many articles; the field dependence 
of the magnetization is presented at the several fixed values of temperatur e 19 ! 20 ' 21 , while the particular details of 
its temperature dependence is overlooked. Hida et al. discussed in the sine-Gordon model 22 the winding numbers 
induced by the chemical potentials, which can be interpreted as the magnetization due to the finite magnetic field in 
the Haldane chain. This mapping describes correctly the h ~ h c case. However, it breaks down at higher magnetic 
fields. Konik et aL— derived the susceptibility as a function of h on the basis of the exact solution of the 

nonlinear sigma (NLcr) model (see M = J x(h)dh in Fig. 5 of Ref 23). Unfortunately, at finite magnetization, the NLcr 
model does not give accurate descriptions as we will discuss shortly. In this way, the comprehensive understanding of 
magnetization in the entire gapless regime h c < h < h s was lacking, which we clarify in the present paper. Various 
universal features are found which also apply to other one-dimensional gapped spin systems. 

We present our findings in the following manner; In Sec. [Til we give the QMC results for the magnetization of 
the S — 1 Haldane chain. In Sec. IHIi we elucidate the universal features of the temperature dependence of the 
magnetization, and explain the numerical results. Our findings on purely one-dimensional systems are also compared 
with the magnon Bosc-Einstcin condensation in three dimensions. 

II. NUMERICAL RESULTS 

Let us first remind of the temperature dependence of the magnetization of the S = 1/2 Heisenberg chain as a 
reference. At present precise numerical data are easily obtained on the basis of the Bethe Ansatz solution (using the 
integrability of the Quantum Transfer Matrh*2i) as shown in fig.[TJfor the various magnetic fields. For any value of the 
magnetic field, as a function of the temperature, the magnetization monotonically increases until it hits a maximum, 
and then monotonically decreases down to zero in the infinite temperature limit. This behavior is qualitatively the 
same as in the classical spin result^. The XXZ chain in the gapless regime also shows a similar behavior. The 
decrease at low temperature is intuitively interpreted as the growth of the antiferromagnetic short range order. 

Next, we present the QMC results for the S = 1 Haldane chain. The QMC simulation is performed at N — 512 with 
the maximum 1200000 Monte-Carlo step using a stochastic series expansio n 10 ' 11 ' 12 code of the ALPS librar y 26 ' 27 ' 28 . 
The convergence in the thermodynamic limit after the finite size scaling is confirmed by the comparison of our results 
with the high temperature expansion results. Figures [2] shows the magnetization in the Haldane gap phase h < h c 
of eq. (p}. As well known, it monotonically decreases and vanishes toward T — > following T _1 exp(— A/T). In 
contrast, the magnetization at h > h c has a characteristic structure. Figures [3]Ja) and[3](b) show the magnetizations 
at h c < h < h s — 4 J (h s is the saturation field). As in the 5 = 1/2 case, the qualitative behavior is consistent with 
that of the classical spin system except at low temperature. However, unlike the S — 1/2 case the magnetization 
minimum is observed at T = T m . The temperature T m , which gives the magnetization minimum, decreases toward 
T = as h is lowered to h c . 
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FIG. 1: The temperature dependence of the magnetizations of the S =1/2 Heisenberg model. The exact classical spin solution 
and the Bethe ansatz results for several choices of (h/J = 0.01 ~ 1.5) are depicted. 
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FIG. 2: The temperature dependence of the magnetization of the S = 1 Heisenberg model for h(< h c ): h = 0.3 J. In this 
regime, the temperature dependence monotonically decreases; T _1 exp( — A/T), especially for T< A. The error bar of each 
point is much smaller than the symbol size. 
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FIG. 3: The temperature dependence of the magnetization of the 5* = 1 Heisenberg model in the gapless phase(/i > h c ) below 
the saturation &eld(h a — 4.0 J) at several values of h/J, (a)0.45-1.5, and (b)2.0-3.8. The exact susceptibility of the classical 
spin Hamiltonian is shown together in (a). The error bar of each point is much smaller than the symbol size. 



III. DISCUSSION 



We discuss the characteristic behavior found in Fig. [3] by use of generic effective theories of the gapped spin chains 
with short-range interactions. Note that the logic given in this section are not restricted to the Haldane chain. 



A. Effective theory near h c and h a 



Excited states above the gap from the S z — ground state generally consists of a triplet massive boson state 
with S z = ±1,0 which usually have repulsive short-range interactions. Low-energy states of these particles may be 
approximated by a non-relativistic dispersion relation: 



E(k) w A + k 2 /2m-hS z . 



(2) 
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where m is a band curvature. If an underlying relativistic theory exists, which is not necessarily the case, vq = yA/m 
corresponds to the relativistic "speed of light" . As a result, the S z = 1 magnon branch intersect with the ground state 
at h = h c = A, a quantum phase transition point. In the gapless regime a quasi long range order appears^, where 
one may integrate out the S z — 0, — 1 magnons which are the higher energy states. Then, in the low density limit of 
the magnons, h — h c , the remaining S z — 1 magnons can be exactly mapped onto the free fermion theory with the 
dispersion eq. ^=21SSL. The low but finite magnon density regime, h > h c , can be still approximately described by 
the free fermion theory since residual interactions are proportional to h — h c . Here, the magnetization is equal to the 
number of the particle as 

T-^r* DW/(e_ '"' <3) 

where e = k 2 /2m, = h — A, and /(e — /i) = l/{exp(/3(e — /i)) + 1} is the Fermi distribution function. The density 
of states follow, D(e) oc 1/y/e, constant, and y/e in one, two, and three-dimensions, respectively. Then, from eq. 
the increase of the magnetization at T — > in Fig [3] can be easily understood by the following argument; the Fermi 
distribution function is symmetric with respect to the Fermi point [i, namely, f(e + 5e — /i) = 1 — /(e — 8e — ji). 
Moreover, f(e — fi) for e > 2[i near T = is negligibly small, since /(e — /x) ~ exp(— /3(e — //)). Therefore, if the density 
of states increases for lower energy e, so does eq. (J3]) at T — > 0. From this observation, we can conclude that this 
decrease of the magnetization in fig. [3] stems from the singularity of the density of states in one-dimensional systems. 
The exact integration of eq. gives the explicit expression, 

T = -^-^-^ W 

where Li„[x] = S^i 2 -'/^" i s the polylogarithm function. Actually, the numerical description of eq. ((H) in fig. 0] 
reproduces well the minimum of M found in fig. [3ja). 

The analogous mapping between the dilute boson and the free fermion holds at h < h s as well. In this case, the 
vacuum state is the fully-polarized state and the low-lying excitations consist of S z = magnon branch instead of the 
S z = 1 one. Then, the magnetization shows a maximum at low temperature corresponding to the minimum of the 
magnon density, just the opposite to that at h > h c . The behavior found in fig.EJb) (e.g h = 3.5,3.8) are consistent 
with this argument. 
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FIG. 4: The free fermion result, eq. Q for h — A = 0.3, is depicted. 



B. Effective theory in the low-energy limit for h c < h < h 3 

The free fermion description has thus succeeded in reproducing the minimum/maximum in fig. [3J However, it is 
valid only slightly above h c and below h s . On the other hand, the TL liquid theory should be applicable for the whole 
gapless regime in one-dimensional systems (h c < h < h s in the present case), albeit only in the low energy limit. In 
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the following, we apply the TL liquid theory to estimate the magnetization. From the conformal field theory, the low 
temperature expansion of the free energy per unit volume is given b y 31 ' 32 

f = ^~^ + 0(T% (5) 

where eo is an ground state energy, vp is the excitation velocity of a fixed point theory, and c is the central charge(c = 1 
for TL liquid). The magnetization is given by the derivative of the free energy with respect to the magnetic field, 

T~*> 

The first and the second term give the magnetization of the ground state and the leading finite temperature correction, 
respectively. This equation indicates that whether the magnetization near T = increases or decreases is determined 
by the sign of the gradient of the velocity with respect to the magnetic field, d^Vp. For the S =1/2 Hcisenberg chain, 
we always have d^v < (see fig. 9 in ref. (3J), which is consistent with Fig.[TJ In contrast, the gapless TL liquid regime 
in the present case has Vp = at both endpoints, h = h c and h — h s , so that Vf must have a maximum in between. 
This indeed gives the characteristic behavior observed in Fig. [3] 

Let us demonstrate the validity of eq. ([6]) in the special case of h > h c , where an explicit evaluation of Vf(/i) is 
available. Following the free-fermion dispersion in eq. © we have 



VF = V °\I A () 

The low temperature expansion of the magnetization given by the free-fermion theory is derived from the Sommerfeld 
expansion^ on eq. Q as, 



M V^U^A) - . /A _ ^ + Q{T% (8) 



L ttvq 12i>o V 2 

with Vq = A/m. These two equations lead to 



M _ v /2A(/i - A) 7T dv Frr , 
L ttvq 0>v 2 f dh 



T 2 + <3(T 4 ), (9) 



which is the special case h > h c of the general result, eq. ([6]). The behavior of M/L at h < h s is analogously obtained 
by replacing vf m eq. ((TJ) with v F oc \fh a — h, following the mapping of S z = 1 at h > h c magnons with S z = 
magnons at h < h s . 

To confirm the above argument we obtain vp, by the DMRG calculation on the 5 = 1 Haldanc chain with periodic 
boundary condition. The energy at the fixed magnetization, M/L, is obtained for discrete values of M = — L with 
system size up to L = 120. One can fit the energy values, E(M/L, L), obtained after the extrapolation of truncation 
error to the finite size scaling: 

E(M/L, L) , , . ttv f (M/L) 1 , T , S 
1 I ' =e (M/L)- F{ 6 1 > — + 0(l/L 3 ), (10) 

where e (M/L) and vp{M/L) denote the energy and velocity at finite magnetization, respectively. Figure[SKa) shows 
the magnetization curve which is obtained from the energy as a function of M/L in its the inset. By use of this 
M/ L — h curve, we plot in Fig. E^b) the velocity as a function of the magnetic field. The velocity increases rapidly just 
above h = h c , takes a maximum at around h ~ J, and then decreases toward h = h s , as anticipated in the previous 
discussion. The reason why the gradient of vp shows a sudden change at h ~ J and 3J is still out of our scheme, 
while the shape of vf in the whole gapless regime suggests that it might be divided into three parts according to the 
change of nature in its excitation spectrum. Going back to Fig. [3] we find several cases where the M — T curves are 
almost flat in a fairly wide temperature range, e.g., h/J — 1.5,2.0. Apparently, this is consistent with the fact that 
the velocity vf takes the maximum (dhVF — 0) at a certain field within this range, as shown in Fig. [5jb) . At this 
field, the leading finite temperature correction starts from 0(T 3 ). Even when h is slightly shifted from the maximum, 
vf depends rather weakly on h. Thus for the range of h around the maximum, magnetization looks almost flat in the 
low temperature regime. 
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FIG. 5: (a) Magnetization as a function of magnetic field of the 5 = 1 Haldane chain in the bulk limit. The inset shows the 
M/L dependence of energy per site obtained by DMRG with periodic boundary condition up to L = 120. The energy at each 
L is extrapolated against the truncation error, where the accuracy is guaranteed to less than the order of 10 -6 in the worst 
case. (b)The magnetic field dependence of velocity derived from the data in (a). The broken lines in (a) and (b) are the free 
fermion values taken from the first term of eq.(9) and from eq.(7) with vo/J = 2.49^2, respectively. 



So far, we have not discussed the effects of irrelevant operators on eq. Generally, Lorentz invariant terms never 
renormalize the spin-wave velocity which corresponds to the "speed of light" . Therefore, to examine whether the 
corrections due to irrelevant operators are required in eq. © , perturbations which break Lorentz invariance should be 
explicitly considered. In principle, irrelevant perturbations could give larger corrections than T 2 , e.g. in the case of 
the S = 1/2 Heisenberg chain, the marginally irrelevant operator gives a logarithmic correction to the zero magnetic 
field susceptibility 3 -. Generally, under the compactification of the bosonic field 4>, the vertex operator exp (ia<f)) (a is 
a constant parameter) and (d x n (f)) n (m and n are integers) can be perturbations around the TL liquid fixed point. In 
the present case with M ^ 0, the former one became less relevant in the low temperature regime, T -C M/J. This 
is because it usually takes a form of exp(zMa;/L) x exp(a0), remaining negligible in the low energy limit. Then, the 
(d x 4>) 3 -term becomes the leading perturbation, since the d x <j) and (d x (j)) 2 '-terms are exactly absorbed into Cq and d 2 4> 
is forbidden by the parity symmetry (x — > — x and <fi — > —(f)). By a standard dimension analysis, this term again leads 
to a T 2 correction to eq. © ■ This correction should be understood as already included in eq. ©, which thus stands 
valid. 



C. Crossover temperature and estimation of the gap 

Let us now return to the discussion to the characteristic temperature that give the magnetization minimum, T m . 
It can be interpreted as a crossover temperature from the TL liquid to the state with non-relativistic dispersion, 
e oc k 2 , which is indicated by the arrow in fig. [6j At T < /i - /i c , the system can be described by a TL liquid with 
linearized dispersions around the two Fermi points (oc \Jh— h c ). However, with increasing T > h — A, the effect of 
band curvature becomes more important due to thermal excitations. 

Here we propose a new way to estimate a gap from T m . Equation (UJ) takes the minimum under the condition, 
2x = Li„ = i/ 2 [— c x ]/Li„ = _ 1/ / 2 [— c x ], where x — n/T m . Its solution, x = x ~ 0.76238, yields 

T m = x (h-A). (11) 

The finite temperature measurement of magnetization at several h > h c in both experiments and numerical calculations 
thus provide a reasonable estimation of A. Since eq. (|11[) consists only of universal constant we can use it without 
any microscopic informations of the system. As a demonstration, we show in Fig. [7] the comparison of our Monte 
carlo results with eq. (fTTj) . The data asymptotically approaches eq. (fTTj) toward h — > A since the deviation at large h 
is due to the interaction effect of the magnons. Note that all the data at h > h c = A in fig. [7] are located below the 
values given by eq. pip . This behavior is rather universal for generic models as follows; the excitations to the higher 
energy modes are enhanced by the repulsive interactions, so that, even in low temperature regime, the band curvature 
effect in the exact results becomes more important than that in the free case. Consequently, the exact value of the 
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crossover temperature, T m , must be suppressed by the repulsive interaction compared to the estimate from the 
free fcrmion theory. 



D. Comments on previous studies 



In light of our results, we comment on the analysis based on the integral )ility of the NLer model, as proposed in 
Ref. l23l . They conclude that the spin-wave velocity monotonically increases with the increasing magnetic field, which 
is contradictory with the results in the present paper. The reason for this discrepancy is that they introduced the 
finite magnetic field after taking the low-energy limit of the zero-field model, which generally does not give a correct 
effective theory above the critical field. To illustrate this point, for a moment, let us consider the 5=1/2 Heisenberg 
antiferromagnetic chain as a simple example. The zero-field effective theory in the low-energy limit is given by the 
TL liquid, namely the free boson field theory with the Lagrangian density, £0 = ^{d^) 2 , where cj> must be periodic, 
4> ~ 4> + 2nR; R is the compactification radius 37 . The application of the magnetic field after taking the the low-energy 
limit can be represented by the term, Cz = ~~^®x4>- This formula appears to be easily handled by "completing the 

square" as Cq + Cz = ^(d^c/)') 2 — j-, where <j)' = <fi— hx/\/2n. Apparently, there is no renormalization of the spin-wave 
velocity nor of R, which is in clear contradiction to the exact Bethe Ansatz solution^. Such fallacy is attributed to 
taking the low-energy limit first at zero-field, as discussed above. Thus one should note that although the irrelevant 
perturbations on the fixed point theory Co can be usually ignored in the low-energy limit, it cannot be in deriving 
the effective theory at a finite field. 

The magnetization minimum has been often discussed in terms of BEC. Although the three-dimensional BEC 
observed in Ref. [l?] gives similar M — T curves to Fig. [3l it differs from the present case in several points; the 
magnetization shows a singular cusp of minimum at the transition temperature, T c cx (h — /i c ) 2 / 3 , (which is however 
smeared in a actual systems 3 ^), whereas T m cx h — h c in our model just marks the crossover. The transverse 
magnetization is finite, namely the off-diagonal long-range order is present in BEC at T < T c , but is absent in our 
one-dimensional system at any temperature. 
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FIG. 6: The phase diagram for the gapped spin chains. Dashed lines indicate the crossovers. The crossover temperature 
between the TLL regime and the quantum critical regime is T ~ h — h c . 



IV. SUMMARY 



We found the characteristic temperature dependence of magnetization of the 5=1 Haldanc chain in the gapless 
regime at finite magnetic field, h c < h < h s . In the vicinity of the endpoints of the gapless regime, h > h c and h < h s , 
the magnetization shows a minimum and maximum, respectively, at the crossover temperature, T = T m . The TL liquid 
theory yields the leading term of magnetization in the low temperature limit as, M(T, h) — M(T = 0, h) cx T 2 . 
We showed that the spin wave velocity, vp(h), vanishes at both h — h c and h s and shows a maximum in between. 
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FIG. 7: The fitting of the QMC results with eq. From this fitting, the gap A can be estimated. The error-bars come from 

a discreteness of the numerical data with respect to temperature. 



Therefore, M(T, h = fixed) at around T ~ 0(T < T m ) undergoes a change from the decreasing function to the 
increasing function of T as the magnetic field increases from h = h c to h s . 

We emphasize that the present discussions are universal (at least) for axial symmetric generic gapped one- 
dimensional spin systems, provided that the field-induced gapless phase can be described by a single-component 
TL liquid. The detailed structure of M and vp reflecting the characteristics of each system does not violate this 
scheme since we have dealt with only the universal features in the low-energy limit. Actually, the conclusions do not 
depend on a particular choice of the effective theory, such as the sine-Gordon model, the nonlinear sigma model, or 
the delta- function Bose gas. 

We also proposed a simple and reliable scheme (represented by eq.(JTTJ) to estimate the gap, A, from the finite 
temperature magnetization measurement. 

To our knowledge, there is no reported experimental results corresponding to our analysis in the present paper. 
However, it might be possible that indeed there are already some corresponding experiments, which have been in- 
terpreted in different ways. Actually, the upturn of magnetization at T = T m found in many experiments had been 
interpreted as a three-dimensional effect (BEC of magnons). However, we argue that it does occur generally in 
purely one-dimensional gapped spin systems as well. Therefore, a care should be taken in the interpretation of the 
magnetization data. To distinguish these two scenarios, one should examine whether or not the transverse magneti- 
zation exists below T m as well as the presence/absence of the thermodynamic phase transition. In reality the actual 
magnetic systems have finite interchain interactions. However, the effects described in the present paper may be 
observable if the interchain interactions are sufficiently weak. We hope that the present work stimulates further study 
on one-dimensional magnetism. 
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